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Abstract 



Due to sphaleron processes in the high-temperature symmetric phase of the stan- 
dard model the cosmological baryon asymmetry is related to neutrino proper- 
5^ , ties. For hierarchical neutrino masses, with B — L broken at the unification scale 

Agut ~ 10"*^^ GeV, the observed baryon asymmetry ub/s ^ 10"^'^ can be naturally 
explained by the decay of heavy Majorana neutrinos. We illustrate this mechanism 
with two models of neutrino masses, consistent with the solar and atmospheric 
neutrino anomalies, which are based on the two symmetry groups SU{5) x U{1)f 
and SU{3)c x SU{3)l x SU{3)r x U{1)f- We also review related cosmological 
bounds on Majorana neutrino masses and the use of Boltzmann equations. 



Contents 



1 Introduction 



2 Cosmological bounds on neutrino masses 

2.1 Chemical equilibrium 

2.2 Relations between B, L and B — L 

2.3 Constraint on Major ana neutrino masses . . . . 

2.4 Primordial asymmetries 

2.5 Supersymmetry 

2.6 R-parity violating interactions 

2.7 Finite mass etl:'ects| 

3 Baryogenesis and neutrino masses 

3.1 Baryogenesis through Majorana neutrino decays 

3.2 Neutrino masses and mixings 

3.2.1 SU{b)xU{l)F 

3.2.2 SU{'i), X SU{'i)L X SU{'i)R x U[1)f . . 

4 Boltzmann equations and scattering processes 

4.1 Boltzmann equations 

O Reaction densities 



5 Baryon asymmetry 

5.1 Delayed decay 

5.2 LeptogenesisI 

5.2.1 ^'(7(5) X (7(1) J. 

5.2.2 SU[^)c^ SU[^)l^ SU[^)rx UJTY^ . 

5.3 Supersymmetric leptogenesis 

Mass scales of leptogenesis and B — L breaking 



B Outlook 



1 Introduction 



The generation of a cosmological baryon asymmetry can be understood as a conse- 
quence of baryon number violation, C and CP violation, and a deviation from thermal 
equilibrium All these conditions can be naturally satisfied in the context of unified 
extensions of the standard model of strong and electroweak interactions. In particular 
the deviation from thermal equilibrium is realized in the out-of-equilibrium decay of 
heavy particles whose mass is related to the mass scale of unification 0. The presently 
observed matter-antimatter asymmetry, the ratio of the baryon density to the entropy 
density of the universe. 
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is then explained as a consequence of the spectrum and interactions of elementary parti- 
cles, together with the cosmological evolution. A general overview of different scenarios 
for baryogenesis can be found in [p|. 
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Figure 1: One of the 12-fermion processes which are in thermal equilibrium in the high- 
temperature phase of the standard model. 

A crucial ingredient of baryogenesis is the connection between baryon number [B) 
and lepton number (L) in the high-temperature, symmetric phase of the standard model. 
Due to the chiral nature of the weak interactions B and L are not conserved |Q . At zero 
temperature this has no observable effect due to the smallness of the weak coupling. 
However, as the temperature approaches the critical temperature of the electroweak 
phase transition B and L violating processes come into thermal equilibrium P]. 



3 



The rate of these processes is related to the free energy of sphaleron-type field config- 
urations which carry topological charge. In the standard model they lead to an effective 
interaction of all left-handed fermions 0] (cf. fig. |l|), 

Ob+l = n {QnqmLilLi) , (2) 

i 

which violates baryon and lepton number by three units, 

Afi = AL = 3 . (3) 

The sphaleron transition rate in the symmetric high-temperature phase has been eval- 
uated by combining an analytical resummation with numerical lattice techniques 0. 
The result is, in accord with previous estimates, that B and L violating processes are in 
thermal equilibrium for temperatures in the range 

Tew ~ 100 GeV <T < Tsph ~ 10^^ GeV . (4) 



Sphaleron processes have a profound effect on the generation of the cosmological 
baryon asymmetry. Eq. ^ suggests that any B + L asymmetry generated before the elec- 
troweak phase transition, i.e., at temperatures T > Tew^ will be washed out. However, 
since only left-handed fields couple to sphalerons, a non-zero value of B + L can per- 
sist in the high-temperature, symmetric phase if there exists a non-vanishing B — L 
asymmetry. An analysis of the chemical potentials of all particle species in the high- 
temperature phase yields the following relation between the baryon asymmetry Yb and 
the corresponding L and B — L asymmetries and Yb-l, respectively [Q, 

Yb = C Yb-l = Ye . (5) 

Here C is a number 0{1). In the standard model with three generations and one Higgs 
doublet one has C = 28/79. 

An important ingredient in the theory of baryogenesis is also the nature of the elec- 
troweak transition from the high-temperature symmetric phase to the low-temperature 
Higgs phase. A first-order phase transition yields a departure from thermal equilibrium. 
Since in the standard model baryon number, C and CP are not conserved, it is con- 
ceivable that the cosmological baryon asymmetry has been generated at the electroweak 
phase transition ^j. This possibility has stimulated a large theoretical activity during 
the past years to determine the phase diagram of the electroweak theory. 

Electroweak baryogenesis requires that the baryon asymmetry, which is generated 
during the phase transition, is not erased by sphaleron processes afterwards. This leads 
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Figure 2: Jump of the Higgs vacuum expectation value at the critical temperature. Com- 
parison of 4d simulations (triangles, squares) with 3d simulations (stars) and perturba- 
tion theory. 



to a condition on the jump of the Higgs vacuum expectation value v = y(f)^(f) at the 
critical temperature |Q, 

— ^ > 1 • (6) 

The strength of the electroweak transition has been studied by numerical and analytical 
methods as function of the Higgs boson mass. The result for the SU{2) gauge-Higgs 
model is shown in fig. ^ [0, where the results of 4d lattice simulations |Tl[, 3d lattice 
simulations [T^ and perturbation theory [jl3[ are compared. For Higgs masses above 
50 GeV this is a good approximation for the full standard model. Since the present 
lower bound from LEP on the Higgs mass has reached almost 110 GeV, it is obvious 
that the electroweak transition in the standard model is too weak for baryogenesis. For 
supersymmetric extensions of the standard model a sufficiently strong first-order phase 
transition can still be achieved for special choices of parameters |1^ . 

For large Higgs masses the nature of the electroweak transition is dominated by non- 
perturbative effects of the SU{2) gauge theory at high temperatures. At a critical Higgs 
mass = 0{my/) an intriguing phenomenon occurs: the first-order phase transition 
turns into a smooth crossover [|15|, |16|, [l^, as expected on general grounds |jlO[- This 
is shown in fig. || [|l^ where X/g"^ ~ mfj / {8m^^) is plotted as function of fi'^/g'^ — 
{3ml,+ml)/{QAV2GFm^){T^-T^)/T\ with = y2m|^/(G^(3m^+m|^)). For small 
Higgs masses one has a first-order phase transition at a critical temperature Tc(m//). At 
the endpoint of the line of first-order transitions, which is reached for mn = the 
phase transition is of second order [IS]. 
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Figure 3: Critical line of first- order phase transitions. 

The value of the critical Higgs mass can be estimated by comparing the W-boson 
mass mw in the Higgs phase with the magnetic mass msM in the symmetric phase. For 
^SM = Cg^T one obtains ||18| , 



/ 3 \i/2 

= j ^ 74 GeV , (7) 



where we have used C ~ 0.35 pnf. Numerical lattice simulations have determined the 



precise value = 72.1 ± 1.4 GeV The estimate (|^) can also be extended to the 
supersymmetric extensions of the standard model, where one obtains for the critical 
Higgs mass < 130 ... 150 GeV m . 



The detailed studies of the electroweak phase transition have shown, that for Higgs 
masses above the present LEP bound of about 110 GeV the cosmological baryon asym- 
metry did not change during this transition, except possibly for a small parameter range 
in the supersymmetric standard model. In particular, the electroweak transition may 
have been just a smooth crossover, without any deviation from thermal equilibrium. In 
this case it's sole effect has been to switch off the B — L changing sphaleron processes 
adiabatically. 

Based on the relation (|^) between baryon and lepton number we then conclude that 
B ~ L violation is needed to explain the cosmological baryon asymmetry if baryogenesis 
took place before the electroweak transition, i.e. at temperatures T > Tew ~ 100 GeV. In 
the standard model, as well as its supersymmetric version and its unified extensions based 
on the gauge group SU{b), B — L is a conserved quantity. Hence, no baryon asymmetry 
can be generated dynamically in these models and one has to consider extensions with 
lepton number violation. 



k 




' 



^0 



Figure 4: Effective lepton number violating interaction. 



The remnant of lepton number violation at low energies is the appearance of an 
effective AL = 2 interaction between lepton and Higgs fields (cf. fig. 



Such an interaction arises in particular from the exchange of heavy Majorana neutrinos. 
In the Higgs phase of the standard model, where the Higgs field acquires a vacuum 
expectation value, it gives rise to Majorana masses of the light neutrinos z/g, z/^ and Ur- 

Lepton number violation appears to be necessary to understand the cosmological 
baryon asymmetry. However, as we shall see in the following section, the required lepton 
number violation can only be weak, otherwise any baryon asymmetry would be washed 
out. These two conditions are the basis of the interesting contraints which the existence 
of the matter-antimatter asymmetry imposes on neutrino physics and on extensions of 
the standard model in general. 

2 Cosmological bounds on neutrino masses 

In the standard model (SM) neutrinos are massless. However, the exchange of heavy 
particles can give rise to an effective lepton Higgs interaction which, after electroweak 
symmetry breaking, generates neutrino masses. The lagrangian describing all fermion- 
Higgs couplings then reads 

Cy = —h^ijdRiQLjHi — h^^jU^iqLjH2 — heije^ikjHi + -fij l^iH2 C kjH2 + h.c. . (9) 

Here q^, um, dm, ki, ^Ri, i = I ■ ■ ■ N, are N generations of quark and lepton fields. Hi 
and H2 are Higgs fields with vacuum expectation values Vi = {H^) 7^ 0. h^, h^, and / 
are NxN complex matrices. For the further discussion it is convenient to choose a basis 



where hu and he are diagonal and real. In the case f 1 ~ f 2 ~ f = y f 1 + f f the smallest 
Yukawa coupling is heii = mg/vi followed by hun = mu/v2- 



J^AL=2 = T^fij lliV C kjV + h.c. . 
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The mixing of the different quark generations is given by the Kobayashi-Maskawa 
matrix V^, which is defined by 

(10) 



Vl 



Here is the diagonal real down quark mass matrix, and the weak eigenstates of 
the right-handed d-quarks have been chosen to be identical to the mass eigenstates. 
Correspondingly, the mixing matrix in the leptonic charged current is determined by 



where 



nil, 
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(12) 



is the diagonal and real mass matrix of the light Majorana neutrinos. 



2.1 Chemical equilibrium 

In a weakly coupled plasma with temperature T and volume V one can assign a chemical 
potential /i to each of the quark, lepton and Higgs fields. In the SM with one Higgs 
doublet, i.e., H2 = Hi = (f, and N generations one has 5N+1 chemical potentials. The 
corresponding partition function is |^| 



Z{fi, T, V) = Tre-^(^-S« ^^^'^ . (13) 

Here f3 = 1/T, H is the Hamilton operator and Qi are the charge operators for the 
corresponding fields. The asymmetry in the particle and antiparticle number densities is 
then given by the derivative of the thermodynamic potential, 

n.-n. = -^^^^^, n{fi,T) = ~lnZ{fi,T,V). (14) 

For a non-interacting gas of massless particles one has 

gT^ \ Pia + O UfSfiif) , fermions , 

Hi — rin = < / C!\ 15 

6 I 2/5/i, + O ((/?/ii)') , bosons . ^ ' 

The following analysis will be based on these relations for (3^i <^ 1. However, one should 
keep in mind that the plasma of the early universe is very different from a weakly coupled 
relativistic gas due to the presence of unscreened non-abelian gauge interactions. Hence, 
non-perturbative effects may be important in some cases. 
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In the high-temperature plasma quarks, leptons and Higgs bosons interact via 
Yukawa and gauge couphngs and, in addition, via the non-perturbative sphaleron pro- 
cesses. In thermal equilibrium all these processes yield constraints between the various 
chemical potentials. The effective interaction (§) induced by the SU{2) electroweak in- 
stantons yields the constraint 

E (3/^</i + f^ii) = • (16) 



One also has to take the SU{3) QCD instanton processes into account |^ which generate 
the effective interaction 

i 

between left-handed and right-handed quarks. The corresponding relation between the 
chemical potentials reads 

J2 (2/Ugi - l^ui - fJ'di) = . (18) 

i 

A third condition, which is valid at all temperatures, arises from the requirement that 
the total hypercharge of the plasma vanishes. From eq. ([T5|) and the known hypercharges 
one obtains ^ 

(j^gi + - fJ'di - f^li - fJ'ei + = • (19) 

The Yukawa interactions, supplemented by gauge interactions, yield relations be- 
tween the chemical potentials of left-handed and right-handed fermions, 

f^qi — fJ-ip — fJ-dj = , fiqi + fiip — fiuj = , fin — fiy^ — fiej = . (20) 

Furthermore, the AL = 2 interaction in (|^) implies 

fJ'ii + fJ'^ = . (21) 

The above relations between chemical potentials hold if the corresponding interac- 
tions are in thermal equilibrium. In the temperature range Tew ~ 100 GeV < T < 
TspH ~ 10^^ GeV, which is of interest for baryogenesis, this is the case for all gauge 
interactions. It is not always true, however, for Yukawa interactions. The rate of a scat- 
tering process between left- and right-handed fermions, Higgs boson and W-boson, 

i^L^ ^ ^rW , (22) 

is r ~ aX'^T, with a = g'^/{4:7i). This rate has to be compared with the Hubble rate, 
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The equilibrium condition T(T) > H{T) is satisfied for sufficiently small temperatures, 

T <Tx^ A^IO^^ GeV . (24) 

Hence, one obtains for the decoupling temperatures of right-handed electrons, up- 
quarks,..., 

Te ~ 10^ GeV , T„ ~ 10^ GeV , . . . . (25) 

At a temperature T ~ 10^° GeV, which is characteristic of leptogenesis, cr = cri, 
Hr = eR2, dR = (Iri, sr = dR2 and ur = uri are out of equilibrium. 

Note that the QCD sphaleron constraint (|18D is automatically satisfied if the quark 
Yukawa interactions are in equilibrium (cf. (pO|)). If the Yukawa interaction of one of the 
right-handed quarks is too weak, the sphaleron constraint still establishes full chemical 
equilibrium. 

Using eq. (|15D also the baryon number density tir = BT'^/6 and the lepton number 
densities = LT'^/Q can be expressed in terms of the chemical potentials. The baryon 
asymmetry B and the lepton asymmetries Li read 

B = J2 C^l^li + l^^i + l^di) , (26) 

i 

Li = 2/i/i + ^lei , L = ^Li . (27) 

i 

2.2 Relations between 5, L and B — L 

Knowing which particle species are in thermal equilibrium one can derive relations be- 
tween different asymmetries. Consider first the most familiar case where all Yukawa 
interactions are in equilibrium and the AL = 2 lepton-Higgs interaction is out of equi- 
librium. In this case the asymmetries Li — B/N are conserved. The Yukawa interactions 
establish equilibrium between the different generations, 

fin = fii , figi = fig , etc. (28) 

Together with the sphaleron process and the hypercharge constraint they allow to express 
all chemical potentials, and therefore all asymmetries, in terms of a single chemical 
potential which may be chosen to be /i;. The result reads 

_2A^ + 3 _6A^ + 1 _2N -1 

1 AN 

= -3/^/ , /^V = ^ 3 /^/ ■ (29) 

The corresponding baryon and lepton asymmetries are 

AN ^ 1AN^ + 9N , , 

-D = III , L = — ui , (30) 

3 Pi , QN + 3 ^ ^ ^ 
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which yields the well-known connection between the B and B — L asymmetries ^1 



B^M±±(B-L). (31) 

Note, that this relation only holds for temperatures T ^ v.\n general, the ratio B /{B — 
L) is a function of v/T [|^, 

Another instructive example is the case where the AL = 2 interactions are in equi- 
librium but the right-handed electrons are not. Depending on the neutrino masses and 
mixings, this could be the case for temperatures above Te ~ 10'* GeV |Q. Right-handed 
electron number would then be conserved, and Yukawa and gauge interactions would 
relate all asymmetries to the asymmetry of right-handed electrons. The various chemical 
potentials are given by = /iei, yUg = Aie2 = • • • = IJ-eN), 

_ 3 _ _ 1 _ 1 

~ ~ lOiV^' ' ^'^ ~ " lOiV^' ' ~ 5iV^' ' 
- 3 _ 1 _ 3 

^' " ~2m^'' ^^"20iV^^' ^^"20iV^^- 
The corresponding baryon and lepton asymmetries are 

1 4N + 3 

B = -f.e, L = ^^^^, (33) 

which yields for the relation between B and B — L, 

2N 

B = (B-L). (34) 

Note that although sphaleron processes and AL = 2 processes are in equilibrium, the 
asymmetries in B, L and B — L do not vanish! 



2.3 Constraint on Majorana neutrino masses 

The two examples illustrate the connection between lepton number and baryon number 
induced by sphaleron processes. They also show how this connection depends on other 
processes in the high-temperature plasma. To have one quark-Higgs or lepton-Higgs 
interaction out of equilibrium is sufficient in order to have non- vanishing B, L and B — L. 
If all interactions in (^ are in equilibrium, eqs. (|2l|) and (|29|) together imply fii = and 
therefore 

B = L = B - L = , (35) 

which is inconsistent with the existence of a matter-antimatter asymmetry. Since the 
equilibrium conditions of the various interactions are temperature dependent, and the 
AL = 2 interaction is related to neutrino masses and mixings, one obtains important 
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constraints on neutrino properties from the existence of the cosmological baryon asym- 
metry. 



The AL = 2 processes described by (M) take place with the rate 



rAL=2(T) = i^^ E <■ (36) 

Requiring T^l=2{T) < H(T) then yields an upper bound on Majorana neutrino masses, 

i:<<(o.2eV (^)"y . (37) 

For typical leptogenesis temperatures T ~ 10^° GeV this bound is comparable to the 
upper bound on the electron neutrino mass obtained from neutrinoless double beta decay. 
Note, that the bound also applies to the r-neutrino mass. However, if one uses for T the 
decoupling temperature of right-handed electrons, Tg ~ IC^ GeV, the much weaker bound 
< 2 keV is obtained . 



Clearly, what temperature one has to use in eq. (|37|) depends on the thermal history 
of the early universe. Some information is needed on what kind of asymmetries may 
have been generated as the temperature decreased. This, together with the temperature 
dependence of the lepton-Higgs interactions then yields constraints on neutrino masses. 



2.4 Primordial asymmetries 

The possible generation of asymmetries can be systematically studied by listing all the 
higher-dimensional SU{3) x SU{2) x f/(l) operators which may be generated by the 
exchange of heavy particles. The dynamics of the heavy particles may then generate an 
asymmetry in the quantum numbers carried by the massless fields which appear in the 
operator. 

For d=5, there is a unique operator, which has already been discussed above, 

{lLiV){lLj¥^) . (38) 

It is generated in particular by the exchange of heavy Majorana neutrinos whose coupling 
to the massless fields is 

huijT^ikj'f ■ (39) 

The out-of-equilibrium decays of the heavy neutrinos can generate a lepton asymmetry, 
which is the well-known mechanism of leptogenesis. The decays yield asymmetries Li — 
B/N which are conserved in the subsequent evolution. The initial asymmetry in right- 
handed electrons is zero. In order to satisfy the out-of-equilibrium condition it is very 
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important that at least some Yukawa couplings are small and that the right-handed 
neutrinos carry no quantum numbers with respect to unbroken gauge symmetries. 

In order to study possible asymmetries of right-handed electrons one has to consider 
operators containing cr. A simple example, with d=6, reads 

{qLilLjWRkem) ■ (40) 

It can be generated by leptoquark exchange (x ~ (3*, 1, 1/3)), 

{>^i,qlCh, + \lul,CeR,)x . (41) 

Note, that (p and x form a 5*-plet of SU{^). In principle, out-of-equilibrium decays 
of leptoquarks may generate a cr asymmetry. One may worry, however, whether the 
branching ratio into final states containing cr is sufficiently large. Furthermore, it appears 
very difficult to satisfy the out-of-equilibrium condition since leptoquarks carry colour. 
Maybe, all these problems can be overcome by making use of coherent oscillations of 
scalar fields |^ or by special particle production mechanisms after inflation. However, 
we are not aware of a consistent scenario for the generation of a Cij- asymmetry. Hence, 
it appears appropriate to take the bound eq. (|37|) on Majorana neutrino masses as a 
guideline and to examine its validity in each particular model. 



2.5 Super symmetry 

At temperatures T where some of the interactions are out of thermal equilibrium the 
effective theory acquires a larger symmetry. Examples discussed above are the AL = 2 
interactions and the Yukawa couplings of the right-handed electron. In the latter case 
the symmetry of the effective theory is generated by QeR-, the number operator of right- 
handed electrons. As a consequence, the chemical potential /ig, i-e. the Ci^-asymmetry, is 
a free parameter whose value influences B and L. 

In the supersymmetric standard model (SSM) one has Yukawa interactions given by 
the superpotential 

W = hd^JD^QjH^ + K,,U^Q,H2 + KijE^L.H, + ^f,,L,H2LjH2 , (42) 

and mass parameters, the supersymmetric mass term ^HiH2 and soft supersymmetry 
breaking scalar masses and gaugino masses, 

Lm = ^mggigL + ^m^WLWL + ]^'rnJ)LhL + h.c. , (43) 

for bino, wino and gluino, respectively. Naturally, the gaugino masses and the /i- 
parameter are of order the gravitino mass, mj ~ /i ~ tjIq. For unbroken supersymmetry, 
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the SSM lagrangian has two chiral U(l) symmetries in addition to those of the SM, an 
axial Peccei-Quinn symmetry and R-invariance. 

Due to the additional fermions, the higgsinos /ili,/il2 and the gauginos, the effec- 
tive instanton induced interactions are modified. For a Weyl fermion ip in the adjoint 
representation of a SU {N) gauge theory the 't Hooft determinant reads, 

Oadj = V^Li • • • ^L2iv • (44) 

For a theory with Weyl fermions in the fundamental and the adjoint representation the 
't Hooft determinant is the product of the two single contributions. Hence, in the SSM 
eqs. (0) and ( ]T7| ) are replaced by 

Ob+L = n {(iLiqiiqiilLi) hLlhL2WLWLWLWL (45) 
i 

and 

0^ = 11 i(lLiqLiU%id%-) gL9L9L9L9L9L , (46) 

i 

respectively. The corresponding relations for the chemical potentials read, 

i^f^Qi + f^li)+ f^hl + f^h2 + = , (47) 

i 

and 

{2figi - fiui - fidi) + GfJ'g = . (48) 

i 

Due to the different Higgs content of the SSM also the zero hypercharge constraint is 
modified, 

i 

+2 E {^^Qi + 2/i„-i - /iJi - l^u - /^el) + -^^^^h2 -f^hi) ■ (49) 

i 

Note, that scalars contribute twice as much as fermions. Furthermore, one obtains from 
the gaugino interactions, 

f^Qi = f^b + = l^qi = /^S + l^qi , (50) 

and therefore, 

^^b = ^^w = ■ (51) 

The gaugino interactions also determine all chemical potentials of scalars in terms of the 
chemical potentials of the corresponding fermions. 
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If all mass terms, i.e., all effects of supersymmetry breaking, are in equilibrium, one 

has 

/^/il + /^/l2 = ' /ife = /id, = /ig = . (52) 

The number of parameters is then the same as in the SM and one obtains essentially the 
same relations between B, L and B — L. 

However, at sufficiently high temperatures, T > Tsb, supersymmetry breaking effects 
will not be in equilibrium. As an example, consider the case that all Yukawa interactions 
are in equihbrium, which implies fin = fii, fiqi = fig, etc. The above relations then 



determine all chemical potentials in terms of ni and jjg [0. The baryon and lepton 
asymmetries 

B = N{2fXg + fiu + iXd) + 2N(2fiQ + fi^ + fi^) , (53) 
L = Ar(2/ii + /ie) + 2Ar(2/i^ + /ig) , (54) 

are then given by, 

lOA^ - 24 

B = -4iV/i; + fig, (55) 

UN'^ + 9N 4A^2^18iV + 3 , , 

^ = 2iV + l 2iV + l 

Clearly, the simple proportionality between B, L and B — L is lost. In particular, one 
obtains for B — L = 0, 

16N'^ + 212A^2 _ 234Ar _ 216 
^ " 22iV2 + 13iV ' 

which yields B = 6 fig for = 3. Hence, at temperatures T > Tsb also an asymmetry in 
gauginos affects the baryon asymmetry. 



2.6 R-parity violating interactions 

Like for Majorana neutrinos, cosmological bounds can also be derived on the strength 



of new baryon and lepton number changing interactions 0, Of particular in- 
terest are R-parity violating interactions in supersymmetric theories which allow single 
production of supersymmetric particles at colliders. 

Consider, as an example, the AL = 1 Yukawa couplings 

W^R = Xi.kDtQjLj, , (58) 

which may appear natural since the Higgs multiplet Hi and the lepton multiplets Lk 
have the same gauge quantum numbers. Requiring Tx < Tew ~ 100 GeV and neglecting 
the masses of scalar quarks and leptons, one obtains from eq. (^) 

A < 10"^ , (59) 
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which would make the interaction (^8]) unobservable in colhder experiments. A detailed 
study P5[ of 2 ^ 1 and 2^2 processes taking finite mass effects into account yields 



bounds slightly less stringent than (|59|). 

It is important, however, that the bound (^) does not apply to all couplings Xijk 
P5| . It is sufficient that B — 3Lj is effectively conserved for a single lepton ffavour. If 
the couplings Xijk have such a ffavour structure and ii a B — 3Lj asymmetry is indeed 
generated a primordial baryon asymmetry can be preserved. Since the corresponding 
couplings A are smaller than all other Yukawa interactions, the generation of a i? — 3Lj 
asymmetry is indeed conceivable. 



2.7 Finite mass effects 



The relation ([15[) between asymmetries and chemical potentials holds for an ideal rel- 
ativistic gas. It is modiffed due to interactions. To leading order in the couplings their 
effect can be expressed in terms of thermal masses [BR p6l B^. In the symmetric, high- 



temperature phase one finds [37 



Hi 



gT^ ( 



3 m^{T) 



J^2 



(60) 



For leptons, for instance, the thermal masses are given by 



3 

32' 

1 2 1,2 

ogi + o^ei 



(61) 



where gi and g2 are the U{1) and SU{2) gauge couplings, respectively. 

The thermal masses also affect the relation between baryon and lepton numbers. 
In particular, one can have B ^ with B — L = if the lepton asymmetry is ffavour 



dependent. In the symmetric phase one finds ^7 

517 



B 



K (Le -Lr + L^~ Lr) 



28447r2 
9 • 10"^ (Le -Lr + L^ 



(62) 



Although this finite mass effect is very small, it can have interesting consequences if 
large ffavour dependent lepton asymmetries are generated in a theory where B — L is 
conserved. 
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3 Baryogenesis and neutrino masses 



3.1 Baryogenesis through Majorana neutrino decays 

As discussed in section 1, baryogenesis before the electroweak transition requires B — L 
violation, and therefore L violation. Lepton number violation is most simply realized 
by adding right-handed Majorana neutrinos to the standard model. Heavy right-handed 
Majorana neutrinos, whose existence is predicted by all extensions of the standard model 
containing i? — L as a local symmetry, can also explain the smallness of the light neutrino 
masses via the see-saw mechanism [^ . 

The most general Lagrangian for couplings and masses of charged leptons and neu- 
trinos reads 

Cy = —h^ijeRilLjHi — h^ijl7RilLjH2 — Ijj^rijl^Ri^RjR + h.c. . (63) 

The vacuum expectation values of the Higgs fields, {Hi) = vi and {H2) = V2 = tan/3 vi, 
generate Dirac masses rrif. and mD for charged leptons and neutrinos, rrif. = h^Vi and 
'^D = h„V2, respectively, which are assumed to be much smaller than the Majorana 
masses M = hr {R). This yields light and heavy neutrino mass eigenstates 

z/ ^ z/^ + , N^un + u%, (64) 

with masses 

m^, ~ -Vjm^— mr.K , tun ^ M . (65) 

Here Vy is the mixing matrix in the leptonic charged current (cf. eqs. (|^) and (pT])). 

The right-handed neutrinos, whose exchange may erase any lepton asymmetry, can 
also generate a lepton asymmetry by means of out-of-equilibrium decays. This lepton 
asymmetry is then partially transformed into a baryon asymmetry by sphaleron processes 



39| . The decay width of the heavy neutrino Ni reads at tree level. 



r,,, = T{N,-^H2 + l) + T (iV, H'^ + r) = ^{hyhlhMi . (66) 

ovr 

A necessary requirement for baryogenesis is the out-of-equilibrium condition F/ji < 



H\t=Mi IMyl' where H is the Hubble parameter at temperature T. From the decay width 



(|6^) one then obtains an upper bound on an effective light neutrino mass p^ , ^ , 

mi = [huhDii-^ cr^Agt,' 77- 

< 10"^ eV. (67) 

Here g^: is the number of relativistic degrees of freedom, Mp = (SnGj^)^^^'^ ~ 2.4 • 
10^^ GeV is the Planck mass, and we have assumed ~ 100, V2 — 174 GeV. More direct 
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Figure 5: Tree level and one-loop diagrams contributing to heavy neutrino decays. 



bounds on the light neutrino masses depend on the structure of the Dirac neutrino mass 
matrix. The bound (|67|) imphes that the heavy neutrinos are not able to follow the rapid 
change of the equilibrium particle distribution once the temperature of the universe drops 
below the mass Mi. Hence, the deviation from thermal equilibrium consists in a too large 
number densities of heavy neutrinos, as compared to the equilibrium density. 

Eventually, however, the neutrinos will decay, and a lepton asymmetry is generated 
due to the CP asymmetry which comes about through interference between the tree-level 
amplitude and the one-loop diagrams shown in fig. In a basis, where the right-handed 
neutrino mass matrix M = (R) is diagonal, one obtains 



r(iVi 
1 



8^ Kht 



11 



i=2,3 
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Ml 



+ g 



'mT 

Ml 



here / is the contribution from the one-loop vertex correction. 



fix) 



1 + x) In 



X 



(68) 



(69) 



and g denotes the contribution from the one-loop self energy [H H, H, which can 
be reliably calculated in perturbation theory for sufficiently large mass splittings, i.e., 

\Mi-Mi\ > ir^-rii, 

gix) = . (70) 



X 



For Ml < M2, M3 one obtains 

3 

'8^ 



£^1 



h,yht 



11 



i=2,3 



Ml 
M 



(71) 



In the case of mass differences of order the decay widths one expects an enhancement 



from the self-energy contribution , although the influence of the thermal bath on this 
effect is presently unclear. 

The CP asymmetry ( p^ ) then leads to a lepton asymmetry pn |, 

Ul - n-r El 



Yr 



K 



9* 



(72) 
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Here the factor k < 1 represents the effect of washout processes. In order to determine 
K one has to solve the full Boltzmann equations In the examples discussed in 

section ^ one obtains k 10^^ . . . 10~^. 

In a complete analysis one also has to consider washout processes. Particularly im- 
portant are AL = 2 lepton Higgs scatterings mediated by heavy neutrinos (cf. fig. ^) since 
cancellations between on-shell contributions to these scatterings and contributions from 
neutrino decays and inverse decays ensure that no asymmetry is generated in thermal 
equilibrium |50[] . 

Further, due to the large top-quark Yukawa coupling one has to take into account 
neutrino top-quark scatterings mediated by Higgs bosons |47, 48]. As we will see in 
section O, these processes are of crucial importance for leptogenesis, since they can 
create a thermal population of heavy neutrinos at high temperatures T > Mi. As the CP 
asymmetry can be interpreted as a mean lepton asymmetry produced per neutrino decay, 
the requested baryon asymmetry can only be generated if the neutrinos are sufficiently 
numerous before decaying. 

Various extensions of the standard model have been considered in early studies of 
leptogenesis In particular, it is intriguing that in the simple case of hierarchical 

heavy neutrino masses the observed value of the baryon asymmetry is obtained without 
any fine tuning of parameters if B—L is broken at the unification scale, Aqut ~ 10^^ GeV. 
The corresponding light neutrino masses are very small, i.e., rriiy^ ~ 3 ■ 10^^ eV, as 
preferred by the MSW explanation of the solar neutrino deficit, and m^.^ ~ 0.1 eV |]58| . 
Such small neutrino masses are also consistent with the atmospheric neutrino anomaly 
E^, which implies a small mass m^^g in the case of hierarchical neutrino masses. This fact 
gave rise to a renewed interest in leptogenesis in recent years, and a number of interesting 



models have been suggested |65|-|81 



3.2 Neutrino masses and mixings 

Leptogenesis relates the cosmological baryon asymmetry and neutrino masses and mix- 
ings. The predicted value of the baryon asymmetry depends on the CP asymmetry (|68D 
which is determined by the Dirac and the Majorana neutrino mass matrices. Depending 
on the neutrino mass hierarchy and the size of the mixing angles the CP asymmetry can 
vary over many orders of magnitude. It is therefore important to see whether patterns of 



neutrino masses |]82[ motivated by other considerations are consistent with leptogenesis. 
In the following we shall consider two examples. 

An attractive framework to explain the observed mass hierarchies of quarks and 
charged leptons is the Froggatt-Nielsen mechanism |^ based on a spontaneously bro- 
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ken U(1)f generation symmetry. The Yukawa couplings arise from non-renormalizable 
interactions after a gauge singlet field $ acquires a vacuum expectation value, 

Q i~\~Q j 



9ij 



(73) 



Here gij are couplings 0{1) and Qi are the U(l) charges of the various fermions, with 
Qq, = —1. The interaction scale A is usually chosen to be very large, A > Aqut- In 
the following we shall discuss two different realizations of this idea which are motivated 
by the atmospheric neutrino anomaly |^^. Both scenarios have a large z/^ — Ur mixing 
angle. They differ, however, by the symmetry structure and by the size of the parameter 
e which characterizes the flavour mixing. 



3.2.1 SU{5) X U{1)f 

This symmetry has been considered by a number of authors. Particularly interesting is 
the case with a non-parallel family structure where the chiral U{1)f charges are different 
for the 5*-plets and the 10-plets of the same family An example of possible 

charges Qi is given in table 1. 

The assignment of the same charge to the lepton doublets of the second and third 
generation leads to a neutrino mass matrix of the form [Q, |85[| , 

(e^ e e\ , 

e 1 1 
[ell) 



(R) 



(74) 



This structure immediately yields a large z/^ — Ur mixing angle. The phenomenology of 

neutrino oscillations depends on the unspecified coefficients 0{1). The parameter e which 

gives the flavour mixing is chosen to be 

($) 1 

= e ~ — . 

A 17 

The three Yukawa matrices for the leptons have the structure. 



(75) 
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Table 1: Chiral charges of charged and neutral leptons with SU{5) x U{1)f symmetry 
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^R2 


^Rl 


^L3 ^L2 ^Ll 




^^^2 


^Rl 


Qi 





1 
2 


5 
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1 5 
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1 
2 


5 
2 



Table 2: Chiral charges of charged and neutral leptons with SU (3)c x SU (3)^ x SU {?>)r x 
U{1)f symmetry |^ . 



Note, that /ig and have the same, non-symmetric structure. One easily verifies that 
the mass ratios for charged leptons, heavy and light Majorana neutrinos are given by 



rrie '■ : m,- ~ e : e : 1 



mi : m2 : 



Ml : M2 : Ms ~ : : 1 
: 1 : 1 . 



(77) 
(78) 



The masses of the two eigenstates z/2 and z/3 depend on unspecified factors of order one, 
and may differ by an order of magnitude [^, |9^. They can therefore be consistent with 
the mass differences Am^^^^ ~ 4- 10~^ — 1 ■ 10~^ eV^ inferred from the MSW solution 
of the solar neutrino problem and Ami — (5 ■ 10~^ — 6 ■ 10"'^) eV^ associated with 



the atmospheric neutrino deficit For numerical estimates we shall use the average 
of the neutrino masses of the second and third family, m^ = {m^^m^^Y^'^ ~ 10~^ eV. 

The choice of the charges in table 1 corresponds to large Yukawa couplings of the 
third generation. For the mass of the heaviest Majorana neutrino one finds 



M, 



^ ~ 10^^ GeV 



(79) 



Since hrss and the gauge coupling of U{1)b-l are C(l), this implies that B — L is broken 
at the unification scale Aqut- 

3.2.2 SU{3)c X SU{3)l x SU{3)r x U{1)f 

This symmetry arises in unified theories based on the gauge group Eq. The leptons e^, 
II and z/^ are contained in a single (1, 3, 3) representation. Hence, all leptons of the same 
generation have the same U{1)f charge and all leptonic Yukawa matrices are symmetric. 
Masses and mixings of quarks and charged leptons can be successfully described by 
using the charges given in table 2 . Clearly, the three Yukawa matrices have the same 
structure^. 



e 

5/2 



e 
e 

ei/2 



^5/2 ^^ 

el/2 

1 



y g5/2 gl/2 



e5/2 ^^ 
el/2 
1 



(80) 



iNote, that with respect to ref. e and e have been interchanged. 
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but the expansion parameter in h,^ may be different from the one in he and hr. From the 
quark masses, which also contain e and e, one infers e ^ P7[ . 

From eq. (|80| ) one obtains for the masses of charged leptons, hght and heavy Majorana 
neutrinos, 

Ml : M2 : M3 ~ : e : 1 , (81) 

mi : m2 : ma ~ e^^ : : 1 . (82) 

Like in the example with SU{5) x U{l)p symmetry, the mass of the heaviest Majorana 
neutrino, 



.,2 

M3 ^ ~ 10^^ GeV , (83) 

ms 



implies that B — L is broken at the unification scale Aqut- 

The Ufj^ — Ur mixing angle is related to the mixing of the charged leptons of the second 



and third generation |87 



sin ~ + e . (84) 

This requires large flavour mixing. 



A / '2 



5) 



In view of the unknown coefficients 0{\) the corresponding mixing angle sinO^,- ~ 0.7 
is consistent with the interpretation of the atmospheric neutrino anomaly as — v-r 
oscillation. 

It is very instructive to compare the two scenarios of lepton masses and mixings 
described above. In the first case, the large — mixing angle follows from a non- 
parallel flavour symmetry. The parameter e, which characterizes the flavour mixing, is 
small. In the second case, the large — mixing angle is a consequence of the large 
flavour mixing e. The V{\)p charges of all leptons are the same, i.e., one has a parallel 
family structure. Also the mass hierarchies, given in terms of e, are rather different. This 
illustrates that the separation into a flavour mixing parameter e and coefficients OiX) is 
far from unique. It is therefore important to study other observables which depend on the 
lepton mass matrices. This includes lepton-flavour changing processes and, in particular, 
also the cosmological baryon asymmetry. 



22 



4 Boltzmann equations and scattering processes 



4.1 Boltzmann equations 

A full quantum mechanical description of baryogenesis has to take into account the in- 
terplay of all processes in the plasma, i.e. decays, inverse decays and scattering processes 
together with the time evolution of the system. Such a quantum mechanical treatment 



may be based either on the time evolution of the density matrix ||^ or on the Kadanoff- 
Baym equations for the Green functions of the system. In the latter case a systematic 
perturbative expansion has recently been obtained which starts from a set of Boltzmann 



equations for distribution functions EM. So far, however, all detailed studies of baryo- 



genesis are based on the Boltzmann equations for number densities. In this section we 
therefore briefly review the basic ingredients of this approach. In particular it is assumed 
that between scatterings the particles move freely in the gravitational field of the ex- 
panding universe and that the interactions are described by quantum field theory at zero 
temperature. 

Consider first a point-particle with mass m > moving freely in a gravitational field. 
Its coordinates in phase space x^, p^^ obey the geodesic equations of motion |9! 

dpf" 



dr 

dx'^ 

li7 



+ T^,yp^ = 0, (86) 
p^ . (87) 



If p^ is the 4-momentum of the particle, the affine parameter r is uniquely determined, 
except for its origin. For a massive particle, m > 0, s = mr is the proper time. With the 
Liouville operator 

the geodesic equations can be written as 

dp'' 



dr 
dx'^ 



L K] , (89) 
. (90) 



dr 

Consider now a gas of non-interacting particles, i.e. of particles whose motion is deter 



mined by the equations (89) and (9C). Due to the absence of interactions, the phase space 
density /^(x,p) has to be constant, i.e. its total derivative with respect to the parameter 
r has to vanish, 

^^^ = 0. (91) 
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If we restrict the Liouville operator (|88D to the mass-shell of these particles, 

where i runs over spatial indices 1 to 3, the coUisionless Boltzmann equation ( PT| ) reads 



95 



Lm[Ui^,p)] = 0. (93) 

In a spatially homogeneous and isotropic universe, described by the Robert son- Walker 
metric, the phase space density can only be a function of t and \p^\- The Boltzmann 
equation is then given by 

L„.[U] = E^^-H\p^fl^ = 0, (94) 



where = ^p^ + m| and H is the Hubble parameter. The equilibrium phase space 
density f^q = exp {—E^/T) solves this equation only in the limit — >• or — oo, 
i.e., only extremely relativistic or non-relativistic particles can be in thermal equilibrium 
in a Roberston- Walker universe. This is due to the fact that the Robert son- Walker metric 



yields 



has no timelike spatially constant Killing vector ||96|. Other solutions of ( plf ) are the Bose- 
Einstein and the Fermi-Dirac distributions of massless particles. 

However, one can come close to thermal equilibrium by including interactions, which 
are described by a collision term C [f^] in the Boltzmann equation: 

E^^-H\P^fl§;^ = C[U] . (95) 
Integrating over the phase space element 

„^^3Hn,^^,l'^ClU], (97) 

where we have made use of the mass shell condition and the spatial isotropy, and is 
the number density. 

The collision term counts the number of collisions a particle i/j undergoes in a time 
and volume element. For the process ip + a + b + . . . —>■ i+j + . . . it is given by [p7| , ^ ^ : 

= J dp^dpadpb . . . dpidpj . . . (27r)^5*^ {p^ + pa + Pb + ■ ■ ■ - Pi - Pj - ■ ■ ■) 

x\M{tP + a + b + ...-^z+j + ...)f UfJb . . . (1 ± /.) (1 ± /,) . . . , (98) 
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where the squared matrix element \A4 {ip + a + b + . . . ^ i + j + . . .)f has to be summed 
over internal degrees of freedom of incoming and outgoing particles. A symmetry factor 
1/n! has to be included if there are n identical incoming or outgoing particles. 

In a dilute gas the factors (1 ± /j), where the upper (lower) sign refers to bosons 
(fermions), can be neglected, and we only have to consider decays ip ^ i + j + . . ., two 
particle scatterings ip + a^i+j + ..., and the corresponding back reactions. Then the 
Boltzmann equation for reads: 

fi^ + 3Hn^ = - [-f{^/j ^ i+ j + ...)- -f{i+ j + ... ^ ^p)] 

id,--- 

- E b{i^ + a^i + j + ...)--f{i + j + ...^tp + a)] . (99) 

a,i,j,... 

Let us consider the decay term first. From the definition (|98D it follows that 
E bi-ip + J + ■■■)- + J + ■■■ ^^)] 



E J dpj . . . (27r)^5^ {p^ - pi 



X 



\M{i,^i + 3 + ...)\^U-\M{i + 3 + ...^n^hf,... 



(100) 



In thermal equilibrium we have fifj ... = because of energy conservation. Hence, 
due to the unitarity of the ^-matrix, this decay term vanishes in thermal equilibrium. 
The same holds true for scattering processes, even if quantum corrections are included 



49, 97, 94 



Further, one can distinguish elastic and inelastic scatterings. Elastic scatterings only 
affect the phase space densities of the particles but not the number densities, whereas 
inelastic scatterings do change the number densities. If elastic scatterings do occur at a 
higher rate than inelastic scatterings one can assume kinetic equilibrium, i.e., the phase 
space density is 



The Boltzmann equation (|99D then takes the form 



(101) 



(102) 



E 
E 

a,i,j,.. 



niUj . . . 

eq eq 



{l+J + ...^ij) 



n^Ua 



^ + J + 



■) eq eq 



ip + a) 



where 7^'' is the space time density of a given process in thermal equilibrium. Note, that 
elastic scatterings no longer contribute to the evolution of n^. 
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It is convenient to replace the particle density by a quantity which is not affected 
by the expansion of the universe. Consider the number of particles in a comoving volume 
element, i.e., the ratio of particle density and entropy density s, 



s 



If the universe expands isentropically. 



(103) 



(104) 



Hence, Y^p is not affected by the expansion, it can only be changed by interactions. 
Further, it is useful to transform to the dimensionless evolution variable z = m^/T. In 



a radiation dominated universe the Boltzmann equation ( |102| ) then finally reads 
sH (m^) dY^ _ 

YY. 



z dz 
'Y 



(105) 



E 
E 



1.1 



Y^Yg 

la 



■i^ J ■ ■ ■ 
yeqyeq 
J-i 1 j • • • 



YHip + a^i + j + ...) 



YY 



Yeqyeq 
J-i ij 



tp + a) 



where H{m^) is the Hubble parameter at T = m^. 



4.2 Reaction densities 

The reaction densities in thermal equilibrium j^*^ can be further simplified. The decay 
width of a particle with energy E is 

r = -r,, , (106) 

where T^s denotes the decay width in the rest system. The reaction density for the decay 
is then given by 

l'S = n^^,^rs: (107) 

here the ratio of Bessel functions is the thermal average of the time dilatation factor 
m/E in eq. 



The decay rate Td for the particle is defined as the number of decays per time 

eg. 



element, i.e., it is given by the ratio of 7^' and the number density n'^- 



D 



K,(z] 



Trs . (108) 



26 



In the low and high temperature hmits the decay rate is given by 

frs , T < 

i.e., due to time dilatation, decays are suppressed at high temperatures. 

In the Boltzmann equation ( |102| ) the expansion term 3Hn^ is responsible for devia- 
tions from thermal equilibrium, whereas interaction terms try to bring the system into 
equilibrium. Hence, the ratio of reaction rates F to the Hubble parameter if is a measure 
for the effectiveness of the interactions. If the reaction rates are too small the system 
will be driven out of equilibrium by the expansion of the universe. 

For inverse decays we know from energy conservation that fi'^fj'^ . . . = . If we 
further neglect possible CP violating effects we have 

M{i) ^ i ^ 3 ^ . . ) = M{% ^ 3 + ■ ■ ■ ^ ^) , (110) 

i.e., in thermal equilibrium, the reaction densities for decays and inverse decays are equal. 
If the decay products are massless, the reaction rate for inverse decays is [^ : 

eg 

r7D = ^r^, (111) 

where = {g/7[^)T^ is the number density of a massless particle species with g internal 
degrees of freedom in thermal equilibrium. In the limit of low and high temperatures we 
now have 



ID 



(112) 



For a two particle scattering process one obtains 



oo / . — 

+ + J d,a(.)y^KW^) , (113) 

where a{s) is the reduced cross section 

2A (s,ml,ml] 
s 

Here, a{s) is the usual cross section as a function of the squared centre of mass energy 
s and A is a flux factor 

X{x', y\ z') = ^[x^-(y + zy] [x^-^y-zy] . (US) 
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If \v\ denotes the relative velocity of the two particles, 



— , (116) 



Ea 

the reaction density for the scattering can also be written as 

^e,^^ + a-.z+j + ...)=n':;!n:^ {a \v\) . (117) 
In the following we shall consider processes with up to three particles in the final state. 

5 Baryon asymmetry 

We are now in a position to evaluate the baryon asymmetry for the patterns of neutrino 



masses discussed in section |3.2| . Before doing that, however, it is useful to consider some 
general aspects of leptogenesis, which are independent of any particular parametrization 
chosen for the neutrino mass matrices. This supplements the qualitative discussion of 



section 3.1. 



5.1 Delayed decay 

A characteristic feature of this baryogenesis scenario is that the generated asymmetry 
depends mostly on the mass parameter 

^1 = (^-^Dii • (118) 

This can be seen in fig. ^ where we have plotted the generated lepton asymmetry as a 
function of mi for three different heavy neutrino masses Mi = 10* GeV, 10^° GeV and 
10^^ GeV. Further, we have assumed a fixed CP asymmetry ei = —10^^ and a fixed mass 
hierarchy for right-handed neutrinos, M| = 10^ Mf, M| = 10^ Mf. 

First one sees that leptogenesis is only possible in a rather narrow range of mi, and 
that the washout processes, parametrized by the factor n in eq. (|72D, are important over 
the whole parameter range, since k < 0.1. If mi is too low, the Yukawa interactions are 
too weak to produce a sufficient number of neutrinos at high temperatures, whereas for 
large mi the washout processes are too strong and destroy any generated asymmetry. 

The asymmetry depends almost only on mi for small mi ^ 10"^ eV, since then the 
generated asymmetry is determined by the number of neutrinos produced at high tem- 
peratures, i.e., on the strength of the processes in which a right-handed neutrino can be 
generated or annihilated. The dominant reactions are decays, inverse decays and scat- 
terings with a top. In the Boltzmann equation ( |106D the reaction densities for heavy 
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Figure 6: Generated B — L asymmetry as a function of mi for Mi = 10^ GeV (dotted 
line), Ml = 10^° GeV (solid line) and Mi = 10^^ GeV (dashed line). The shaded area 
shows the measured value for the asymmetry. 

neutrino production, 'Jprod, all give contributions proportional to mi at high tempera- 
tures T > Ml, 

For large mi ^ 10~^ eV on the other hand, the neutrinos reach thermal equilibrium at 
high temperatures, i.e., the generated asymmetry depends mostly on the influence of 
the lepton number violating scatterings 7al=2 at temperatures T^Mi. In contrast to 
eq. ( |119| ) the lepton number violating processes mediated by a heavy neutrino behave at 
low temperatures like 

with ^ 

= (^-^l^)^-^-^ ' (121) 

where we have neglected interference terms. One therefore expects that the generated 
asymmetry becomes smaller for growing neutrino mass Mi and this is exactly what one 
observes in fig. |^. 
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Eq. (|120|) can also explain the small dependence of the asymmetry on the heavy 
neutrino mass Mi for mi ^ 10"'^ eV. The inverse decay processes which take part in pro- 
ducing the neutrinos at high temperatures are CP violating, i.e., they generate a lepton 
asymmetry at high temperatures. Due to the interplay of inverse decay processes and 
lepton number violating 2 — > 2 scatterings this asymmetry has a different sign compared 
to the one generated in neutrino decays at low temperatures, i.e., the asymmetries will 
partially cancel each other. This cancellation can only be avoided if the asymmetry gen- 
erated at high temperatures is washed out before the neutrinos decay. At temperatures 
T > Mj the lepton number violating scatterings behave like 

^0'..-ocM.i:f S|. (122) 

Hence, for fixed heavy neutrino mass hierarchy, the wash-out processes are more efficient 
for larger neutrino masses, i.e., the final asymmetry should grow with the neutrino mass 
Ml. The finally generated asymmetry is not affected by the stronger wash-out processes, 
since for small rhi the neutrinos decay late, where one can neglect the lepton number 
violating scatterings. 



5.2 Leptogenesis 

We can now evaluate the baryon asymmetry for the two patterns of neutrino mass 
matrices discussed in sections p.2.1| and p.2.2| . Since for the Yukawa couplings only the 
powers in e are known, we will also obtain the CP asymmetries and the corresponding 
baryon asymmetries to leading order in e, i.e., up to unknown factors 0{1). Note, that 
for models with a U{1)f generation symmetry the baryon asymmetry is 'quantized', i.e., 
changing the U{1)f charges will change the baryon asymmetry by powers of e |[88[| . 



5.2.1 SUib) X U{1)f 

In this case one obtains from eqs. ( |68D and (^) 



From eq. ([72D, ~ 1/300 ([75|) and ~ 100 one then obtains the baryon asymmetry, 

Yb ~ K 10"^ . (124) 

For K ~ 0.1... 0.01 this is indeed the correct order of magnitude! The baryogenesis 
temperature is given by the mass of the lightest of the heavy Majorana neutrinos, 

Tb ~ Ml ~ e^Mg ~ 10^° GeV . (125) 
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Figure 7: Time evolution of the neutrino number density and the lepton asymmetry in the 
case of the SU{5) x U{1)f symmetry. The solid line shows the solution of the Boltzmann 
equation for the right-handed neutrinos, while the corresponding equilibrium distribution 
is represented by the dashed line. The absolute value of the lepton asymmetry is given 
by the dotted line and the hatched area shows the lepton asymmetry corresponding to the 
observed baryon asymmetry. 



This set of parameters, where the CP asymmetry is given in terms of the mass 
hierarchy of the heavy neutrinos, has been studied in detail |5^. The generated baryon 
asymmetry does not depend on the flavour mixing of the light neutrinos, in particular 
the — Vj. mixing angle. The solution of the full Boltzmann equations is shown in fig. |^ 
for the non-supersymmetric case. The initial condition at a temperature T ~ lOMi is 
chosen to be a state without heavy neutrinos. The Yukawa interactions are sufficient to 
bring the heavy neutrinos into thermal equilibrium. At temperatures T ~ Mi this is 
followed by the usual out-of-equilibrium decays which lead to a non-vanishing baryon 
asymmetry. The final asymmetry agrees with the estimate ( |124D for k ~ 0.1. The dip in 
fig. 7 is due to a change of sign in the lepton asymmetry at T ~ Mi, as discussed in the 
previous section. 



5.2.2 5f/(3)e X SU{?,)l X SU{?,)r x U{1)f 



In this case the neutrino Yukawa couplings ([801 ) yield the CP asymmetry 

^1 ~ t|- , (126) 
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Figure 8: Solution of the Boltzmann equations in the case of the SU{3)c x SU{3)l x 
SU{3)rX U{1)f symmetry. 



which correspond to the baryon asymmetry (cf. ([72|) ) 

Yb ~ ft 10"^ . (127) 

Due to the large value of e the CP asymmetry is two orders of magnitude larger than in the 
case with SU{5) x ^7(1)^ symmetry. However, washout processes are now also stronger. 
The solution of the Boltzmann equations is shown in fig. The final asymmetry is again 
Yb ~ 10~^ which corresponds to /t ~ 10"'^. The baryogenesis temperature is considerably 
larger than in the first case, 

Tb ~ Ml ~ e^Ma ~ 10^^ GeV . (128) 



The baryon asymmetry is largely determined by the parameter mi defined in eq. ( pTj ) 
48|. In the first example, one has mi ~ fn^. In the second case one finds mi ~ m^. 



Since m^, and are rather similar it is not too surprizing that the generated baryon 
asymmetry is about the same in both cases. 



5.3 Super symmetric leptogenesis 



Without an intermediate scale of symmetry breaking, the unification of gauge couplings 
appears to require low-energy supersymmetry. Therefore, we are going to briefly review 
supersymmetric leptogenesis |]53| , ^ |60[] in the following. As in the non-supersymmetric 
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Figure 9: Decay modes of the right-handed Majorana neutrinos and their scalar partners 
in the supersymmetric scenario. 
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Figure 10: Same as fig. \^for the supersymmetric version of leptogenesis. 



case, a full analysis of the mechanism including all the relevant scattering processes is 
necessary in order to get a reliable relation between the input parameters and the final 
asymmetry. 

The supersymmetric generalization of the lagrangian is the superpotential 



(129) 



where, in the usual notation. Hi, H2, L, E'^, N'^ and R are chiral superfields describing 
spin-0 and spin-| fields. The basis for the lepton fields can be chosen as in the non- 
supersymmetric case. 

The heavy neutrinos and their scalar partners can decay into various final states 
(cf. fig. pp. At tree level, the decay widths read. 



Vrs{N{ ^ 1 + H2) = r,,(iVi ^l + h' 
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The CP asymmetry in each of the decay channels is given by |^ 
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Figure 11: Solution of the Boltzmann equations for supersymmetric leptogenesis for the 
parameters given by the SU{5) x U{1)f symmetry. Y^^ stands for the equilibrium dis- 
tribution of neutrinos Ni, and Fjv^ and Yi^ are the solutions for the neutrino number 
and the sum of scalar neutrino and scalar anti-neutrino number per comoving volume 
element, respectively. The generated lepton asymmetries in leptons and scalar leptons 
are denoted by Yj;^^. and Y^^, whereas Yi_ stands for the difference of scalar neutrino and 
scalar anti-neutrino numbers. 



where 

fix) = iix) [in + ^] . (133) 

It arises through interference of tree level and one-loop diagrams shown in fig. In 
the case of a mass hierarchy, Mj ^ Mj, the CP asymmetry is twice as large as in the 
non-supersymmetric case. 

Like in the non-supersymmetric scenario lepton number violating scatterings medi- 
ated by heavy (s)neutrinos have to be included in a consistent analysis. In the supersym- 
metric case a large number of processes contribute which can easily reduce the gener- 
ated asymmetry by two orders of magnitude. Similarly, the large number of (s)neutrino 
production processes makes leptogenesis possible for values of mi smaller than in the 
non-supersymmetric case ||6^ . 



In fig. we have plotted the generated lepton asymmetry as function of mi for 
three different values of Mi, where we have again assumed the hierarchy M| = 10^ Mf, 
M| = 10^ Mf and the CP asymmetry ei = — 10~^. Fig. |10| demonstrates that in the 
whole parameter range the generated asymmetry is significantly smaller than the value 
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Figure 12: Solution of the supersymmetric Boltzmann equations in the case of the 
SU{3)c X SU{3)l X SU{3)rX U{1)f symmetry. 

4 ■ 10^^ which one obtains, if one neglects lepton number violating scattering processes. 
Baryogenesis is possible in the range 

10"^ eV ^ mi ^ 5 ■ 10"^ eV . (134) 

Comparing non-supersymmetric and supersymmetric leptogenesis one sees that the 
larger CP asymmetry and the additional contributions from the sneutrino decays in the 
supersymmetric scenario are compensated by the wash-out processes which are stronger 
than in the non-supersymmetric case. The final asymmetries are of the same order in 
the non-supersymmetric and in the supersymmetric case. 

Like in the non-supersymmetric scenario it is interesting to see whether a parameter 
choice which gives light neutrino masses and mixing angles compatible with experimen- 
tal results, can also explain the baryon asymmetry. To this end we again considered the 
two scenarios discussed in sections ^.2.1| and |3.2.2| . The time evolution of the (s)neutrino 



number densities and the (s)lepton asymmetries in the case of the SU{5) x U{1)f sym- 



metry is shown in fig. |TT[ Although the CP-asymmetry is now larger than eq. (|123|) by a 
factor of two, the wash-out processes are more effective than in the non-supersymmetric 
case, and the final asymmetry is again of order 10~^, corresponding to a washout factor 
K ~ 0.03. Similarly, in the case of the SU{3)c x SU{3)l x SU{3)r x U{l)p symmetry 
(cf. fig. [T^) the final asymmetry is again of order 10~^. The baryogenesis temperature 
Tb ~ 10^° GeV in the SU{5) x U{1)f model is consistent with the constraint from the 
allowed gravitino abundance [p8| . 
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Leptogenesis can also be considered in extended models which contain heavy SU{2)- 
triplet Higgs fields in addition to right-handed neutrinos Decays of the heavy 

scalar bosons can in principle also contribute to the baryon asymmetry. However, since 
these Higgs particles carry gauge quantum numbers they are strongly coupled to the 
plasma and it is difficult to satisfy the out-of-equilibrium condition. The resulting large 
baryogenesis temperature is in conffict with the 'gravitino constraint' | |99|| . 



5.4 Mass scales of leptogenesis and B — L breaking 

In the previous section we have considered baryogenesis in two models with hierarchical 
neutrino masses but rather different ffavour structure. In both cases the Majorana masses 
of the right-handed neutrinos turned out to be rather large, i.e., M, > 10^° GeV, which 
implies large mass scales for the leptogenesis temperature and B — L breaking. 

In order to understand the model dependence of this result it is instructive to study 
the case of two generations for which bounds on light and heavy neutrino masses can be 
easily obtained (cf. PH ). 



Without loss of generality the heavy Majorana neutrino mass matrix can always be 
assumed to be diagonal with positive real eigenvalues, 



M = . (135) 




Similarly, in terms of the matrix mD = im£,Vy (cf. (^)) the light Majorana neutrino 
mass matrix reads 

= ml^rn, = ly"^^ ^ J , (136) 
where the eigenvalues rrii are again real. With 



= (1 + Tj^) . (138) 




one obtains 



Note that a . . .d are complex. Since mi and 1112 are real there is only one independent 
phase. A complete set of parameters for the two neutrino mass matrices is given by 
mi, m2. Ml, M2 and rf = pe*". 

The out-of-equilibrium condition for the decay of the heavy neutrino implies an 
interesting upper bound on the mass of the lightest neutrino m_ = min{mi,m2}. From 
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< H\t=Mi one obtains for the effective neutrino mass mi, 



With 



one is led to 



mi = ^ ^ <mo = 10-3 eV . (139) 



~ mi + m2|r/r /./in^ 
"^1 = — TT"^ — T\ — ' (140) 
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< 10^3 eV . (141) 

Note, that this bound is only accurate up to a factor 0(1), as the values for rhi show 
which we have obtained for the two models discussed above. This difference illustrates 
the accuracy to which the out-of-equilibrium condition holds. 

The upper bound (|141D can always be satisfied by making the neutrino Yukawa cou- 
plings very small. However, leptogenesis also requires a sufficient amount of CP asym- 
metry. From (^) one obtains (Mi < M2), 

II 3 1 , , _t I ^1 

" I^(^;^''"("^"-)-'m2 

(142) 

We now assume a maximal phase, i.e., rf = ip, a = 7r/2, and |£:i| > eo ~ 10~®, which is 
the smallest value required by the observed baryon asymmetry Yb ~ 10^^°. This yields 
for the heavy neutrino mass Mi, 

M ^ ^2 (l + p^)^/^(mi + m2p) 

3 — mfl p 

Since the right-hand side has a minimum at p = (mi/m2)^/'^, one obtains a lower bound 
on Ml, 

/ 2/3 , 2/3\3/2 

Ml > ^ . (144) 

For quasi-degenerate neutrinos the bound (|141|) implies mi ~ m2 < mo ~ 10^'^ eV. This 
yields 

Ml > 10^^ GeV eo ~ 10^° GeV , (145) 

independent of our knowledge about the solar and atmospheric neutrino anomalies. For 
hierarchical neutrinos all masses have to satisfy the upper bound on the electron neutrino 
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mass obtained from /^-decays since the mass differences infered from tlie athmospheric 
and solar neutrino anomalies are much smaller. Hence, one has m_ ^ m_|_ < leV which 
implies 

Ml > 10^^ GeV eo ~ 10^ GeV . (146) 

The bounds ( |145| ) and ( p.46|) confirm the expectation, in accord with the explicit models 
discussed above, that the smallness of the light neutrino masses suggested by experimen- 
tal data as well as the out-of-equilibrium condition requires heavy Majorana neutrinos 
far above the electroweak scale. The bounds (|145|) and (|146D can be relaxed in the case 
I Ml — M2I = 0{ri). Then a resonant enhancement of the CP asymmetry can occur |^|. 
In this case the lower bound Mi > 10^ GeV has been obtained in a recent analysis ||100 



In the context of a grand unified theory the large Majorana masses will be generated 
by a Higgs mechanism which breaks B — L spontaneously. The corresponding vector 
bosons have to be sufficiently heavy such that the out-of-equilibrium condition for the 



decaying Majorana neutrino A^i is not violated by processes like NiNi Z' ^ ee*^ p7 |. 
Assuming Higgs couplings hr < 1 and a gauge coupling 0{1) one obtains the lower bound 
on the scale oi B — L breaking 

A._.~M..>10"(^) . (147) 

In the case of hierarchical heavy Majorana neutrinos the scale oi B — L breaking is 
naturally identified with the grand unification scale, as the models illustrate which we 
discussed in the previous section. 



6 Outlook 

Detailed studies of the thermodynamics of the electroweak interactions at high tempera- 
tures have shown that in the standard model and most of its extensions the electroweak 
transition is too weak to affect the cosmological baryon asymmetry. Hence, one has to 
search for baryogenesis mechanisms above the Fermi scale. 

Due to sphaleron processes baryon number and lepton number are related in the high- 
temperature symmetric phase of the standard model. As a consequence, the cosmological 
baryon asymmetry is related to neutrino properties. Generically, baryogenesis requires 
lepton number violation, which occurs in many extensions of the standard model with 
right-handed neutrinos and Majorana neutrino masses. In detail the relations between 
i?, L and B — L depend on all other processes taking place in the plasma, and therefore 
also on the temperature. 

Although lepton number violation is needed in order to obtain a baryon asymmetry, 
it must not be too strong since otherwise any baryon and lepton asymmetry would be 
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washed out. Hence, leptogenesis leads to stringent upper and lower bounds on the masses 
of the light and heavy Majorana neutrinos, respectively. 

The solar and atmospheric neutrino deficits can be interpreted as a result of neutrino 
oscillations. For hierarchical neutrinos the corresponding neutrino masses are very small. 
Assuming the see-saw mechanism, this suggests the existence of very heavy right-handed 
neutrinos and a large scale oi B — L breaking. 

It is remarkable that these hints on the nature of lepton number violation fit very 
well together with the idea of leptogenesis. For hierarchical neutrino masses, with B — L 
broken at the unification scale Agut ~ 10^^ GeV, the observed baryon asymmetry 
ub/s 10~^° is naturally explained by the decay of heavy Majorana neutrinos. The 
corresponding baryogenesis temperature is Tb ~ 10^° GeV. 

The cosmological baryon asymmetry provides important constraints on neutrino 
masses and mixings. In addition the connection between lepton fiavour and quark fiavour 
changing processes can be studied in unified theories. In supersymmetric models implica- 
tions for the mass spectrum of superparticles can be derived from the cosmological bound 
on the gravitino number density. It is intriguing that the baryogenesis temperature Tb is 
of the same order as the supersymmetry breaking scale A5. It will be interesting to see 
to what extent further theoretical work and experimental data will be able to identify 
or exclude leptogenesis as the dominant source of the cosmological matter- antimatter 
asymmetry. In order to obtain a theory of leptogenesis one also has to go beyond the 
classical Boltzmann equations and to treat the generation of the baryon asymmetry fully 
quantum mechanically, which presents a challenging problem of non-equilibrium physics. 
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